Spin and interaction effects on charge distribution and currents in one-dimensional 
conductors and rings within the Hartree-Fock approximation 
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Using the self-consistent Hartree-Fock approximation for electrons with spin at zero temperature, 
we study the effect of the electronic interactions on the charge distribution in a one-dimensional 
continuous ring containing a single 8 scatterer. We reestablish that the interaction suppresses the 
decay of the Friedel oscillations. Based on this result, we show that in an infinite one dimensional 
conductor containing a weak scatterer, the current is totally suppressed because of a gap opened 
at the Fermi energy. In a canonical ensemble of continuous rings containing many scatterers, the 
interactions enhance the average and the typical persistent current. 



PACS numbers: 72.I0.Fk, 73.20.Dx 

The effects of electronic interactions on characteristic 
properties, such as charge fluctuations, persistent cur- 
rents (PC's) and the conductance of electronic systems 
are very rich and interesting. jj| They strongly depend 
on the strength and range of the interactions, [||]|||] 011 
the dimensionality of the system, and on whether the 
space is discrete or continuous. ^J,[j) Approximate calcu- 
lations, like Hartree-Fock, introduce a great deal of sim- 
plifications, but at the same time many effects may be 
washed out. However, approximate calculations may be 
used to shed more light on specific problems, while keep- 
ing in mind their limitations. In this work we consider 
e-e interactions within the self-consistent Hartree-Fock 
approximation (SCHFA) for electrons with spin at zero 
temperature. For simplicity we assume an equal number 
of electrons of opposite spin states. 

Our aim is to study numerically the interaction effects 
on the charge distribution and the currents in continuous 
one-dimensional (ID) isolated rings and open conductors 
containing a single 5 scatterer, P,[7|,p|,p|p0| as well as on 
the PC's in rings containing many scatterers. [ pl]p2| Even 
within the Hartree-Fock approximation we recover the 
bosonization Band the density-matrix renormalization- 
group result: We show that for a single scatterer in 
a ring the repulsive electronic interaction suppresses the 
decay of the charge oscillations. Based on this we show, 
as a central result, that for an open conductor with a 
weak scatterer the electronic conduction at the Fermi en- 
ergy vanishes because of Bragg reflection coexisting with 
a gap at the Fermi energy. The zero conduction of the 
interacting system was obtained in Refs. |7|,||,D by exact 
and by renormalization group calculations. Within the 
first iteration of the SCHFA, it was shown |^,^| that an at- 
tempt to explain this result by a scattering perturbation 
series is inadequate because of logarithmic divergences 
of the transmission amplitude at the Fermi energy in all 



orders of the series. 

Although the dissipative conductance of the infinite 
conductor is suppressed by the interactions, the PC in 
a ring is not. This is because the conductance depends 
on the properties of the levels close to the Fermi energy 
but the PC is a thermodynamic property that depends on 
the response of all occupied levels. |llj Moreover, we show 
that once many scatterers are considered, the interactions 
not only do not suppress the PC, but even enhance it. 
We write the HF equation for electrons in a ring of 
radius R with angular coordinate and energy units 
h 2 /m e R 2 = 1 (we drop the background term) as 
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The twisted boundary condition ip(0 + 2tt) = 
tp(9) exp(i27T0/<^o) accounts for a flux tf> threading the 
ring. 4>o = hc/e is the flux quantum. Vdi S (#) is the dis- 
order potential which may include a single or many scat- 
terers. The first (second) integral term is the Hartree 
(Fock) term. The electronic wave functions ^>i(0) = 
ipl(6) exp (—i9(j)/(j)o) in the Fock term are 2tt periodic for 
any value of flux. / enumerates the energy levels together 
with the spin state si . N e is the total number of electrons 
in the ring. The cutoff e 2 allows (as in quasi ID) using 
the 3D Coulomb law S and makes the integrations fi- 
nite. The square of the distance between the particles is 
defined || by (6 - 9') 2 = min[|0 - 6>'| 2 , (2tt -\0- 0'\) 2 ]. 
In Eq. (1), ro = eh 2 /m e e 2 denotes the Bohr radius with 
dielectric constant e (to be distinguished from the cutoff 
e). We define the coefficient g = R/tq to be the inter- 
action strength, g ~ 1 corresponds to semiconductors. 
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Because the sum J2vli (O')Vi'(O) represents al- 
most a closure relation we replace, as discussed in Refs. 
|L4| and [Jl5|], the integrodifferential equation (Eq. (Q)) 
by an ordinary Schrodinger equation that we solve self- 
consistently: 
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Here V c g(9) is given by 
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where Re stands for real part. The spin degree of free- 
dom is very important. For spinless electrons the inter- 
action effect is weak because the Fock and Hartrec terms 
tend to cancel each other due to opposite signs and sim- 
ilar absolute values. Taking into account the spin de- 
gree of freedom, the Hartree term is twice as large as 
the exchange term. Then the former dominates V e ff and 
enhances screening; therefore we expect the interaction 
effects to be stronger for electrons with spin. This ex- 
plains the importance of considering spin |l6[ in order to 
understand disordered interacting systems. 

We begin by studying the interaction effect on the 
charge oscillations in a ring with a single scatterer, 



V dia (0) = XS(9). 
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For a strong scatterer, A > Ef (Ef is the Fermi energy), 
the interaction effect on the decay of the charge oscil- 
lations is weak and may even be neglected because the 
scatterer is dominating. For a weak scatterer, A << Ef, 
at the level of the SCHFA we recover the numerical result 
of Ref. | |To| ] based on the density-matrix renormalization 
group: With increasing repulsive interaction g the decay 
of the Friedel oscillations is suppressed (indicating also 
the reliability of our SCHFA) . Figure |l| depicts the decay 
rate for the strongest interaction for which the SCHFA 
still converges. As Fig. 2 shows, the effective potential 
tends to be periodic with half a Fermi wavelength peri- 
odicity. Both (direct and exchange) terms tend to have 
this periodicity which is independent of the interaction 
strength. This behavior holds for a larger number of elec- 
trons on a ring for a given constant charge density. 

The above results may be used to study the effect on 
the charge oscillations and on the conduction in the case 
of a single weak scatterer XS(x) embedded in an infinite 
ID conductor (x is the spatial coordinate). 

For noninteracting electrons the orthogonal wave func- 
tions, with a given spin state, are MM 
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FIG. 1. Charge oscillations per spin, along half circum- 
ference of a ID continuous ring, induced by a weak sin- 
gle scatterer XS(8) (A = 3.8). Thin (bold) symbols stand 
for noninteracting (interacting) electrons. The interactions, 
at self-consistency of the Hartree-Fock calculations, tend to 
make p(6) periodic and to minimize p(0) further. The inter- 
action strength is g — 3, and the flux threading the ring is 
(j>/4>0 = 0.05 (see text). The total number of electrons per 
spin is 40. The curves are the estimations by the indicated 
formulas. 

with k > and A = X/(h 2 /m e ) having units of inverse 
length. r(k, X) = -iX/(k + iX) and t(k, A) = k/(k + iX). 
Because of time-reversal symmetry and the symmetry 
of the potential under coordinate inversion, t' = t and 



-{r/t)*t = 
Ap(x) = 2 



The fluctuating density per spin is 
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Im takes the imaginary part of the exponential integral 
Ei, z = 2\x\{Kf + iX), and Kf is the Fermi wave vector. 
A,o(0) = — 2A tan -1 (if// A) is a minimum. For KfX > 1, 
the asymptotic expansion of Ei implies 
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For the SCHFA the initial V e g is calculated using the 
wave functions of noninteracting electrons. The charge 
fluctuations define the Hartree potential 



V H (x) = g s 
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FIG. 2. The effective potential along the half circumfer- 
ence of the ring of Fig. 1. The thin (bold) symbol is the 
noninteracting (interacting) result. Note the clear tendency 
of V e ff(9) at self-consistency to become periodic, and to screen 
the scattering potential XS(9) (not shown). 



where g s — 1 (2) for electrons without (with) spin. Our 
approximated Fock potential is 



V F (x) 
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dx' 
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Clearly, V g — V H +V F is a function of \x\ , and will change 
during the iterations until self-consistency is reached. V e ff 
is small due to a weak coupling constant (g ~ 1). At 
this point we invoke an approximate self-consistency by 
adopting a suppression [HJhJ of the decay of the Friedel 
oscillations, as was demonstrated above to be valid in the 
SCHFA. We substitute by hand the limit 0] S = in 



Ap(x) 



| sin(2.K>|a;| + t?/) 



(11) 



for Eqs. (|9|) and (|1C|), assuming that this yields a V c h 
close to that from the SCHFA. To carry out the integra- 
tion [in Eqs. (|) and @], we use a cutoff that allows 
contributions only from \x — x'\ > e. This cutoff is equiv- 
alent to that used in Eq. (0). For KfX ;§> 1, we then 
obtain, up to an additive constant, 

V eff (x) = U[g s sm(2K f x + ri f ) ~sm(AK f x + T] f )}, (12) 

where U = —2^/1^(2^/6), and Cj is the cosine integral. 
Equation ( p^ ) shows that V e s has two periodicities: A//2 
from the direct potential (A/ = 2n/Kf), and A//4 from 
the exchange potential. The overall periodicity is given 
by the larger period. The electrons at the Fermi energy 
exactly obey the Bragg condition [JlTj for total reflection, 
i.e., 




FIG. 3. The interaction effect on the statistics of the 
sample persistent current (in units of the PC of the clean 
ring of noninteracting electrons). Thin (bold) symbols 
stand for noninteracting (interacting) electrons for the same 
canonical ensemble of 201 realizations. The interaction 
(g — 1, 0/0o = 0.325) enhances the persistent current. 
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All the states with \k\ < Kf remain practically unaffected 
by the weak and periodic V e s. Note that consistently 
with Eq. ( p"3| ) there is a gap JlTj of order U at the Fermi 
energy. Thus the current vanishes at the Fermi energy. 

For a ring with a weak scatterer the interaction will 
not destroy the PC even if the current at the Fermi en- 
ergy (assuming a large ring) is totally suppressed by the 
periodic effective potential. This follows from the fact 
that all occupied levels contribute to the PC, except at 
Ef, where Eq. (13) is assumed to be relevant. 

In the following we will consider the general case of a 
large number of scatterers in a ring. Figure 2 already 
shows the importance of screening for a single scatterer. 
This indicates that screening is of particular relevance for 
the case of many random scatterers: 



N, 



(14) 



Here the location and strength of the jth scatterer are 
uniformly distributed in (0, 27r) and (—A, A), respec- 
tively. N s is the total number of scatterers in a ring. 
For the numerics we use A = 14 (in scaled units). 

The characteristic features of disordered noninteract- 
ing samples were, e.g., discussed by Imry and Shiren. |1£] 
For noninteracting electrons, the localization length 14] 
at Ef = 200 is £ ~ ir/2. This should reduce the average 
current in open conductors by a factor ~ 1/50. The av- 
erage sample PC of noninteracting electrons was reduced 
by factor ~ 1 /40 which is slightly greater than predicted 
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FIG. 4. The interaction (g = 1.75) enhances the average 
sample persistent current (in units of the PC of the clean 
ring) and introduces a preffered diamagnetic current direc- 
tion. Thin (bold) symbols represent noninteracting (interact- 
ing) electrons for the same 150 realizations. 



for open conductors. The typical sample PC, -J (I 2 ), was 
reduced by factor ~ 1/10 which indicates the importance 
of a statistical study. The fixed total number of electrons 
in a ring was 32 ± 4. 

Figure 3 shows the interaction effect on the sample 
PC statistics for an interaction coupling constant g = 1. 
The interaction reduces the peak, centered at zero, while 
broadening the distribution. Furthermore, the distribu- 
tion gains more weight at negative values of the PC, 
which indicates a diamagnetic tendency. We found that 
the interaction enhances the typical PC (by factor ~2); 
the average PC is neither enhanced nor suppressed. Fig- 
ure 4 shows that for increasing interaction, g = 1.75, 
the average PC is also enhanced by factor ~2. Figures 
3 and 4 both show a clear tendency of the interaction 
to enhance the PC for electrons with spin. For spin- 
less electrons the PC was found JlJ] to be rather unaf- 
fected by interaction. This shows an essential difference 
between models of electrons with or without spin. In 
addition, a clear difference between tight-binding mod- 
els and continuous model s [|19| becomes apparent: In the 
former it was concluded, pof using exact diagonalization 
and the SCHFA, that switching on the e-e interaction in 
the regime of moderate disorder further suppresses the 
PC because of the Mott transition. J2]j] In continuous 
models this transition appears to be irrelevant, since the 
continuous models correspond to tight-binding models at 
very low fillings. Jl^fl 

In conclusion, using the SCHFA in one-dimension, we 
showed the tendency of the electronic interaction to build 
up nondecaying charge oscillations in a ring containing a 
single weak scatterer. Adopting this result for an infinite 
conductor implies a periodic effective potential. The elec- 



tronic conduction was shown to vanish, because of Bragg 
reflection that coexists with a gap at the Fermi energy. 
This shows that, even in the HF limit the influence of the 
interactions on the Friedel oscillations and on conduction 
in one-dimension, calculated by exact and renormaliza- 
tion methods, may be reproduced. In rings the PC is not 
suppressed by the interaction. It is even enhanced in the 
case of many moderate scatterers due to screening. To 
demonstrate these effects, we considered the spin degree 
of freedom and used continuous conductors and rings. 
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